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Abstract. We present results of a comprehensive mult i- frequency study of the radio galaxy B3 J2330-f 3927. The 
1'.'9 wide radio source, consisting of 3 components, is bracketed by 2 objects in our Keck K—hand image. Optical 
^ . , and near-IR Keck spectroscopy of these two objects yield z = 3.087 ± 0.004. The brightest (K = 18.8) object 

' has a standard type II AGN spectrum, and is the most likely location of the AGN, which implies a one-sided jet 

radio morphology. 

Deep 113 GHz observations with the IRAM Plateau de Bure Interferometer reveal CO J = 4 — 3 emission, which 
peaks at the position of the AGN. The CO line is offset by 500 kms"'^ from the systemic redshift of the AGN, 
but corresponds very closely to the velocity shift of an associated HI absorber seen in Lya. This strongly suggests 
that both originate from the same gas reservoir surrounding the AGN host galaxy. Simultaneous 230 GHz 
interferometer observations find a 3x lower integrated fiux density when compared to single dish 250 GHz 
observations with MAMBO at the IRAM 30m telescope. This can be interpreted as spatially resolved thermal 
dust emission at scales of 0'.'5 to 6". 

Finally, we present a r < 1.3% limit to the HI 21 cm absorption against the radio source, which represents the 
seventh non-detection out of 8 z > 2 radio galaxies observed to date with the WSRT. 

We present mass estimates for the atomic, neutral, and ionized hydrogen, and for the dust, ranging from 
M(HI) = 2 X lO'' M0 derived from the associated HI absorber in Lya up to M(H2) = 7 x lO^'^ M© derived from 
the CO emission. This indicates that the host galaxy is surrounded by a massive reservoir of gas and dust. The 
if— band companion objects may be concentrations within this reservoir, which will eventually merge with the 
central galaxy hosting the AGN. 



Key words. Galaxies: individual: B3 J2330-(-3927 - galaxies: active - galaxies: formation - galaxies: emission lines 
- radio lines: galaxies - cosmology: observations 



1. Introduction 

The study of high redshift galaxies has for a long 
time been limited to the optical/near-IR and radio do- 
mains. Recent technological progress has now opened 
up other wavelength regimes (sub-mm, X-ray and even 
7-ray) to high redshift studies. Together with the dis- 
covery of opt i cally- s electe d 'Lyman break' ga laxies {e.g. 
ISteidel et all Il99(l Il999t lOiichi et all l200lh . this has 
much broadened our knowledge of high redshift objects, 
which was previously based on galaxies containing bright 
AGN. Similarly, our knowledge of AGN has also bene- 
fited from this multi-wavelength approach. The fraction 
of galaxies hosting AGN can now be estima ted from deep 
blank-field X-ray a n d opt ical surveys [e.g. iBareer et al.l 
l200l[ ISteidel et all I2fl02f) . providing an estimate of the 
duration of their active phase. 

The presence of an AGN introduces a number of com- 
plications when studying the host galaxy. This is especially 
true for Type I (broad-lined) AGN, where the host galaxy 
can only be seen after a careful subtraction of the dominat- 
ing AGN emission. This problem does not manifest itself 
in Type II (narrow-lined) AGN, where the bright central 
emission is conveniently shielded by optically thick mate- 
rial. We can therefore use the AGN signatures to pinpoint 
the host galaxies, while still having a relatively unbiased 
view of the stellar population. 

The downside of using type II AGN is that their 
signatures are much fainter at optical wavelengths, and 
thus intrinsically more difficult to detect out to high red- 
shifts. Surveys for distant type II AGN have therefore re- 
lied on manifestations of the AGN at other wavelengths. 
Historically, the most successful technique has been the se- 
lection using powerful radio synchrotron emission from the 
AGN, which has lead to th e discovery of almost 150 radio- 
loud type II AGN at z > 2 l)De Breuck et al.U2000bl). with 
the m ost distant known at z = 5.19 llvan Breugel et al.l 
Il999|) . This predominance of radio selection has lead to the 
almost unique identification of high redshift type II AGN 
as 'high redshift radio galaxies' (HzRGs, z > 2). However, 
recently radio-quie t type II AGN at z > 2 have also been 
found in submm l lvison et al.l 1 19981 but see Vernet fc 
Cimatti, 2001), X-rav jStern et all l200| iNorman et all 
l2002(l and optical surveys l|Steidel et al.Ll2002l) . 

Despite these new developments, radio galaxies remain 
amongst the best studied galaxies at high redshift, because 
they represent some of the most luminous stellar systems 
known out to z > 5 (e.g. iBest. Longair. fc R.ottgerineL 
Il998t IPe Breuck et al. . l2002l) . For example, deep optical 
spectropolarimetry can reveal a wealth of info rmation on 
the s tellar population of the host galaxy {e.g. iDev et all 
Il997j) . the presence of a scattered quasar component {e.g. 



IVernet et alll2001a|) . and the kinematics, ionization mech- 
anism and nietallicity of the ex t ended emission line re- 
gions (e.g. IVillar-Martm et all l2000t iDe Breuck et all 
l2000bl) . Their emission lines (especially Lya) are often 
bright enough to derive informati on on the dens i ty and 
metallicity of associated absorbers llvaii Oiik et all Il997at 
iBinette et all l2000Hjaryis et alllioo^ 

The radio selection has the advantage of being insen- 
sitive to the dust properties, unlike optical selection tech- 
niques, which may miss heavily obscured systems. Such 
objects may be identified by their strong dust emission 
at far-IR w avelengths. Ear l y obse rvati ons at 800 /xm and 
1.2 mm by iDunlop et all |l994|) and llvisonl l)l995|) de- 
tected strong thermal dust emission in the radio galax- 
ies 4C 41.17 {z = 3.8) and 8C 1425-1-635 {z = 4.25), 
suggesting that sever al HzRGs do indeed co ntain signifi- 
cant amounts of dust. I Archibald et aLll)200l[l observed 47 
1 < z < 4. 4 radio galaxies at 8 50 ^m using the SCUBA 
bolometer ([Holland et allll999t) on the JCMT. Even after 
correcting for the strong negative fc— correction effects on 
the Rayleigh- Jeans tail of the thermal dust spectrum {e.g. 
iBlain fc Longaiii Il993() . they found a substantial increase 
in the 850 /xm luminosity with redshift. This is confirmed 
by iReuland et all l)2003cj) . who detect more than half of 
the sources with SCUBA in a sample of 15 radio galax- 
ies at z > 3. Such a redshift dependence is not seen for 
quasars, whose mm and submm luminosit y remains al- 
most unchanged between z ^ 2 and z ~ 4 I Omont et all 
I2OO3I: IPriddev et all l2003(l . From a matched sample of 
1.37 < z < 2.0 quasa r s and radio galaxies selected at 
151 MHz. lWillott et all l(2002t) found that the quasars are 
actually > 2 times brighter at submm wavelengths than 
the radio galaxies. They also find an anti-correlation be- 
tween the far-IR luminosity and radio size, and interpret 
this as a synchronization between the jet-triggering event 
and processes controlling the far-IR luminosity. A possible 
explanation for the difference between the (sub)mm red- 
shift evolution of radio galaxies and quasars could then 
be that HzRGs are more massive, leading to a longer 
(sub)mm-luminous phase due to a larger supply of gas. 

The relatively bright submm fiux densities of HzRGs 
imply that they contain dust masses of 10® — IO^Mq, as- 
suming standard conversion factors (see §4.2). The pres- 
ence of the AGN introduces an ambiguity on the mech- 
anism powering this dust emission: direct heating by 
the AGN o r heating by massiye stars formed in star- 
bursts (e.g. ISanders fc Mirabel! ll99(iHCarim et allbOOOt 



lOmont et ^.'."^20011 IWill'ottretl^l.L^f200^^ ^ Altho"h the 
AGN could easily introduce sufficient energy to produce 
this emission, an increasing number of authors have ar- 
gued that the starbursts arc th e dominating mechanism 
{e.g. iRowan-Robinson. 2000: Ta,dhunter et a,1.. ,20021. A 
strong argument in favour of the starburst origin is the 
detection of dust and CO emission from regions s everal 
tens of kpc away from t he AGN l|Omont et all Il996l: 
IPapadopoulos et alll2000|) . The high gas masses inferred 
from the CO detections could then be the likely reser- 
voirs feeding these massive starbursts. The presence of 
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CO in such reservoirs implies that they must be (at least 
partially) metal emiched, which is consistent with the re- 
cent detection of ionized metal lines at scales of several 
tens of kpc in the z = 2.92 2 radio galaxy MRC 0943-242 
l)Villar-Martm et alllioO^ . These metals may have been 
deposited by previous merger even ts, or by galaxy-scale 
outflows {e.g. iDawson et all 12002(1 . 

These large dust and molecular gas masses are rem- 
iniscent of the often huge L y g ha l oes seen around 
HzRGs (e.g. Ivan Oiik et all 1199(4 IVenemans et all 
l2002t IVillar-Martin et all l200a iReuland et all l2003b() . 
Moreover, if— band observations show that HzRGs 
are amongst the r nost massive galaxies kno wn at 
each redshift (e . g. iBest. Longair. fc Rottgerinel Hi 



iDe Breuck et all l2002l) . However" to date, the relation- 
ship between the different gas and dust components and 
the AGN host galaxy has remained elusive. What is the 
role in the galaxy formation process of the extended Lya 
haloes and the associated absorbers seen in the rest-frame 
UV emission lines? Are they composed of primordial ma- 
terial, or have they been enriched by material from the 
host galaxy? Are the dust and CO reservoirs in any way 
linked to these? 

In this paper, we report the discovery of dust and CO 
emission in the z = 3.087 radio galaxy B3 J2330-I-3927. 
This is only the third z > 3 radio galaxy where CO 
has been detected, after 4 C 60.07 and 6C 1908-^722 
{z = 3.791 and z = 3.532; IPapadoDOiilos et all l200nl) 
and despite several int ensive searches l|Evans et all Il996t 
Ivan Oiik et aliri997lj) . Moreover, the CO emission is co- 
incident in velocity space with associated HI absorption 
in the Lya line. This shows for the first time that these dif- 
ferent components are indeed physically related, and trace 
a massive reservoir surrounding the AGN host galaxy. 

Throughout this paper, we use a cosmology with a 
A— dominated Universe with Hq = 65/i65 km s~^ Mpc~^, 
r^M = 0.3 and 17a = 0.7. At z = 3.087, the luminosity dis- 
tance is Dl = 28.34/155^ Gpc = 8.74 x lO^^h"^ m, and 1" 
corresponds to 8.2 kpc. For comparison, in a flat Universe 
with Ho = 65/i65 km s"^ Mpc'^ = 1-0 and fl^ = 0.0, 
Dl = 19.07/i^5^ Gpc, and 1" corresponds to 5.5 kpc. 



2. Observations and data reduction 

2.1. Source selection 

IVigotti et aD l|l999^ defined the B3-VLA sample with the 
aim of constructing a homogeneous spectral database for 
a sample of 1049 radio sources, containing flux densities 
in the range 151 MHz to 10.6 GHz. From this sample, 
we selected 104 sources with steep optical spectra {a < 
— 1; 5* oc i^"), which remained undetected on the Palomar 
Observatory Sky Survey. Such sources are known to be 
good candidates to find HzRGs {e.g. iDe Breuck et alj 
l2000al) . We first identify the host galaxies of these candi- 
dates in sensitive optical or near-IR images, followed by 
optical spectroscopy to determine their redshift. In the fol- 
lowing, we report on one of these sources, B3 J2330-I-3927. 



2.2. Optical and Near-IR Imaging and Spectroscopy 

On UT 1995 August 31, wc obtained a 600s i?— band im- 
age of B3 J2330-t-39 27 with the Low Resolution Imaging 
Spectrograph (LRIS: l0ke et allll995|) on the Keck I tele- 
scope. Conditions were photometric with I'.'O seeing. We 
reduced the image using the standard techniques in the 
NOAO IRAF package. We determined the astrometric so- 
luti on using 19 st ars also detected in the USNO 2.0 cata- 
log l|Monet et all 1^)98) . The uncertainty in relative opti- 
cal/radio astrometry is dominated by the absolute uncer- 
tainty of the optical reference frame, which is ~0'.'4 (90% 
confidence limit; Deutsch, 1999). Figure ^ presents a find- 
ing chart, based on the Keck i?— band image. 




23 30 2S 27 26 25 24 23 22 

RIGHT ASCENSION (J2000) 



Fig. 1. i?-band image of B3 J2330+3927 obtained with 
Keck/LRIS. The NVSS position of the radio source is 
marked with a cross. The coordinates of the offset star 
are listed in Tabled 

On UT 1996 September 10, we obtained a 1200s spec- 
trum with LRIS on the Keck I telescope. We used the 
300 £ mm^^ grating blazed at 5000A with the 1" wide slit 
resulting in a r esolution of - lOA (FWHM). The slit was 
oriented at PA=168°, and passes through objects a and b 
(see Fig. EJ. We corrected the spectra for overscan bias, 
and flat-fielded them using internal lamps taken after the 
observations. We extracted a one-dimensional spectrum 
using a 1" wide aperture. Because the conditions were not 
photometric, we initially flux calibrated the spectrum with 
an archival sensitivity function. We then determined the 
i?— band magnitude from the spectrum using the IRAF 
task sbands, and adjusted the flux scale to i? = 24.04, 
which is the 1" diameter aperture magnitude from our 
photometric i?— band image. The main uncertainty in this 
procedure stems from the accuracy of the slit positioning 
on the peak of the i?— band emission. We estimate the 
resulting spectrophotometry to be accurate up to ~30%. 



4 



Carlos De Breuck et al.: CO emission and HI absorption in B3 J2330+3927. 



On UT 1997 September 12, w e obtained near-IR imag- 
ing with NIRC l)Matthews et ali ll994) on the Keck I tele- 
scope. Conditions were photometric with 0'.'45 seeing in 
i^5— band. The total integration time of the ifg— image 
was 1140s. Observing procedures, calibration and data 
reduction techniques w ere similar to those described by 
iDe Breuck et al.l ()2002j) . To determine the astrometric so- 
lution, we used 16 stars also detected in our LRIS i?— band 
image. We estimate the relative near-IR/radio astrometric 
uncertainty to be '-^O'.'S. 

Following the i^5— band imaging, we also obtained 
near-IR spectroscopy with NIRC. We used the grl20 grism 
in combination with the i/if— filter. The position angle 
on the sky was 302°, and includes only object a. This 
position angle was dictated by the need for a sufficiently 
bright guide star in the offset guider. The dispersion was 
'^60 A/pix, and the slit width was 0'.'68, resulting in a res- 
olution of --240 A (FWHM). We obtained 10 exposures of 
180s each. We discarded two exposures which were use- 
less due to poor seeing and/or high sky background vari- 
ability. The total integration time is thus 1440s. We re- 
duced the spectra follo wing the procedures described by 
iDe Breuck et all ((l998l) . We extracted a one-dimensional 
spectrum using a 1"2 wide aperture. Because no near-IR 
line-lamps were available, we used the dispersion formula 
listed in the NIRC manual to wavelength calibrate the 
spectrum. This does not allow an independent estimate 
of the uncertainties in the calibration, but a comparison 
of the 3 detected emission lines with the optical spec- 
trum suggests it is accurate up to ~ 70 A. To correct 
for atmospheric absorption lines, we divided the spec- 
trum by the atmospheric transmission curve for Mauna 
Kea, produced using the program IRTRANS4^. To flux 
calibrate the spectrum, we used a spectrum of the AO 
star IID220750, which was obtained immediately subse- 
quent to our B3 J2330-f 3927 observations. We reduced 
this spectrum in the same way as B3 J2330-f3927. To 
construct a sensitivity function, we assumed IID220750 
has a blackbody spectrurn with T^s =9850 K implying 
V ~ K ^ -0.03 ll.IohnsonLlTifil . Using V = 7.01 from 
the SIMBAD database, this provides a first order flux cal- 
ibration. To refine the absolute flux calibration, we deter- 
mined the Xg— band magnitude from the spectrum us- 
ing the IRAF task sbands, and adjusted the flux scale to 
K = 19.18, which is the 1" diameter aperture magnitude 
from our photometric Ks— hand image. 

We corrected both the op tical and near-IR spectra for 
Galac tic reddening using the' CardelH. Clavton. fc Mathid 
l|l989t) extinction curve wi th Ay = 0.39, as determined 
from the dust maps of ISchlegel. Finkbeiner. fc Davij 
(|l99i). 



^ This curve is available on the UKIRT world wide web 
pages. 



2.3. VLA radio imaging 

To determine the morphology of the radio source, we 
observed B3 J2330-H3927 for 1 hour with the VLA 
l|Napier. Thompson fc Ekersl Il983|) on UT 2002 March 
30. We used the A-array at a frequency of 8.46 GHz 
and a bandwidth of 50 MHz, resulting in a resolution 
of 0'.'35 X 0'.'24. Due to technical problems, no primary 
flux calibrator was observed. We used an observation of 
3C 286 from the next day to calibrate the data. However, 
the absolute flux calibration cannot be considered reliable, 
and we therefore only use this observation to determine 
the source morphology. We used the standard data reduc- 
tion techniques in the NRAO AIVS package. Because 
the source is too faint, we did not self-calibrate the data. 
Figure El shows the resulting 8.46 GHz image. 




23 30 25.1 25.0 24.9 24.8 24.7 24.6 

RIGHT ASCENSION (J2000) 



Fig. 2. 8.46 GHz VLA map of B3 J2330-I-3927 overlaid on 
our Keck band image (subject to a 0'.'5 uncertainty 
in the registration). We identify 3 components a, b, and 
c, which may be part of the same physical system (the 2 
small peaks west and north-west of object a are bad pix- 
els). The cross marks the position in the NVSS catalogue, 
which is within 0'.^2 from the p o sition found in the 1.4 GHz 
A-array data of IVigotti et"al1 l)l999|) . The FWHM of the 
synthesized 8.46 GHz VLA beam is shown in the bottom 
left corner. 

2.4. IRAM interferometric Observations 

B3 J2330-f 3927 is one of the brightest 850 fj,m sources 
de tected in the SCUBA survey of z > 3 radio galaxies 
bv lReuland et alJ l)2003c(l . Because gas and dust emission 
are often related, we selected B3 J2330-I-3927 as a prime 
target to search for redshifted CO emission in HzRGs. 

We observed B3 J2330-H3927 with the IRAM Plateau 
de Bure Interferometer during 8 individual sessions be- 
tween 2001 May and 2002 February. We used 3 differ- 
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ent configurations (D, C and B) with 4 to 6 antennae, 
depending on availability. The total on-source observing 
time was ~27 hours (excluding calibrations), with an ex- 
cellent sampling of the uv-plane (Fig. The resulting 
112.798 GHz synthesized beam size is 3'.'2 x 2'.'3 at posi- 
tion angle 56°. The correlator setup for the 3 mm receivers 
involved 4 x 160 MHz modules covering a total bandwidth 
of 560 MHz tuned in single sideband (SSB) mode. 

The central frequency was originally tuned to 
112.834 GHz to observe CO (J = 4-3) {v,^^t = 
461.040 GHz) at z = 3.086, the redshift of the hya 
line. However, after an initial assessment of the data from 
the first session, we re-centered the central frequency to 
112.780 GHz {z = 3.092) to better cover the fuU extent 
of the detected emission line, which appeared to be offset 
from the optical redshift (see below). 

We simultaneously observed at 230.538 GHz (1.3 mm) 
in double side-band (DSB) mode to probe the spatial ex- 
tend of the thermal dust emission. However, in the fol- 
lowing, we consider only those 1.3 mm data that were ob- 
tained when the atmospheric water vapour content was 
< 3 mm. The total 1.3 mm integration time is there- 
fore only 17.7 h. The synthesized beam size at 1.3 mm 
is 1'.'9 X 1'.'5. 

After phase and amplitude calibration, we combined 
the data with different central frequencies into an ex- 
panded data-set of 271 spectral channels of 2.5 MHz each, 
calculating the proper statistical weight of each channel in- 
dividually. The outer channels of the spectral cube there- 
fore have larger uncertainties (see Fig. |HJ • 




V [KMo Wavelsngth] 



Fig. 3. The u - v plane coverage for B3 J2330+3927 
observed at 113 GHz with the IRAM Plateau de Bure 
Interferometer. Note the excellent coverage obtained with 
3 different configurations using 4 to 6 antennae. 



2.5. IRAM Single Dish Bolometer Observations 

To determine the total 1.3 mm flux density, 
B3 J2330-I-3927 was observed in service mode with 
the 117-channel Ma x Planck Bolometer array (MAMBO; 
iKrevsa et all Il998|) at the IRAM 30m telescope on Pico 
Veleta, Spain. The effective bandwidth centre for steep 
thermal spectra is ^^250 GHz. We observed the source 
with the array's central channel, using the standard 
'ON-OFF' observing sequence, in which the secondary 
mirror chops in azimuth by 32" at a rate of 2 Hz. We 
adopted a calibration factor of 35000 counts per Jansky, 
which we estimate is accurate to within 20%. 

The observations were split in 2 sessions: 5x16 sub- 
scans of 12s on UT 2002 October 30 and 8x16 subscans 
on UT 2002 November 2, resulting in a total on-source 
integration time of 21 minute s. We reduce d the data with 
the MOPSI software package llZvlkaLll998l) . including sub- 
traction of the correlated sky noise. This yields a target 
flux of S'250GHz = 4.8 ± 1.2 mJy. 

2.6. WSRT 21 cm line Observations 

We observed B3 J2330-f3927 with the Westerbork 
Synthesis Radio Telescope (WSRT) on several occasions 
between UT 2001 July 1 and 2002 September 10 to search 
for redshifted HI 21 cm (1420.4057 MHz) absorption 
against the synchrotron emission from the radio source. 
Unfortunately, a significant fraction of these observations 
were rendered useless due to external radio frequency in- 
terference (RFI), or by an unstable ionosphere. In the fol- 
lowing, we shall therefore only use the highest quality ob- 
servations, with a total integration time of 22 hours. 

We used the DZB backend with a 10 MHz bandwidth 
centered on 347 MHz and 128 channels giving a velocity 
resolution of ^ 16kms^^ before Hanning smoothing the 
data. 

We made a line cube using uniform weighting, giving a 
resolution of 61" x 36" at PA=0 (north-south). The noise 
in the single channel is ^ 1.7 mJy /Beam/channel (after 
Hanning smoothing) . We have also obtained a continuum 
image of the field. At the observed frequency and resolu- 
tion of the WSRT observations the continuum emission of 
B3 J2330-I-3927 is unresolved with a peak flux of 405 mJy, 
consistent with the 325 MHz WENSS and 365 MHz Texas 
flux densities (see also Fig. . 

3. Results 

3.1. Identification and spectroscopy of the 
optical/near-IR components 

Figure El shows the location of the radio source on the 
Ks— band image. The radio source, consisting of 3 com- 
ponents with a largest separation of 1"9, is located ex- 
actly in between two objects detected in both the R— and 
i^s— band images. We also detect a resolved component 
in the i^g— band image 3"5 southwest of the radio source. 
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Table 1. Astrometry and photometry (in 2" apertures) 
of the optical/near-IR components of B3 J2330+3927. 



Object 


RA(J2000) 


DEC(J2000) 


-R(mag) 




a 


23''30"24?85 


39°27'12''6 


23.08±0.12 


18.79±0.03 


b 


23f> 30^24=94 


39°27'10''3 


24.18±0.22 


20.76±0.11 


c 


23''30™24?70 


39°27'08''9 


>25 


20.70±0.10 


Offset 


23''30'" 23^71 


39°27'22''8 


15.3 


14.4±0.2 


CO peak 


23''30'"24=84 


39°27'12''2 







which may be related. TableHlists the positions and pho- 
tometry of these 3 components, which we label as a, &, 
and c. 

The optical spectroscopic slit included objects a and 
b, while the near-IR slit included only object a. Figures 01 
and El show that object a has the classical spectrum of a 
type II AGN, showing narrow emission lines and a faint 
continuum emission. We only detect a single, weak emis- 
sion line in object b (Fig. 01), which can only be interpreted 
as Lya at z = 3.107; we can exclude the alternative iden- 
tification as [on] A 3727 at z = 0.339, because in this 
case, we would have also detected the [O III] A 5007 line. 
At z = 3.107, the [OIII] A 5007 hne may contribute sig- 
nificantly to the _ftr— band flux, but object b cannot be 
a reflection nebula, as emission is also detected in the 
i?— band, which is free of strong emission lines. Object 6 
therefore appears to be a companion object, offset by 
-1-1500 km s^^with respect to the systemic redshift of ob- 
ject a. 



Figure IHl shows the 1-dimensional optical and near-IR 
spectra of B3 J2330-I-3927 (only object a is included in 
the apertures). We detect 4 emission lines in the opti- 
cal spectrum (see Table EJ , and 3 emission lines in the 
near-IR spectrum. Because the low resolution and uncer- 
tain wavelength calibration of the near-IR spectrum, we 
only use the optical spectrum to determine the systemic 
redshift of B3 J2330-h3927. We adopt the redshift of the 
Hell A 1640 line, z = 3.087 ± 0.004, because this is the 
brighest non-resonant line. 

The HK spectrum allows us to determine the contri- 
bution from emission lines to the iiT— band magnitude of 
object a. We flnd that the [O III] A 5007 line contributes 
36% (0.5 magnitudes) to the Ks— hand flux (the II/3 line 
is much weaker and can be safely ignored). This value 
is comparable to the 45% cont ribution in the z = 3.594 
radio galaxy MG J2144-fl928 (jArmus et al.L Il998l) . and 
to the 34 % average contribution in a sample four ra- 
dio galaxi es where the [O III] A 5 007 line falls in the 
K-hand l|Eales fc Rawlingsl Il993(l . After the emission 
line contribution is removed, the 64 kpc if —band mag- 
nitude {K — 18.79 ± 0.03) is consistent with those fo r 
powerful HzRGs at this redshift l|van Breugel et a'llll998|) . 
This shows that the inclusion of the strong [O III] A 5007 
line can increase the scatter in the K — z di agram, but it 
is un likely to be a dominant factor (see also l.Tarvis et all 
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Fig. 4. Two-dimensional spectrum, centered on Lya. 
Object 6 is a companion object, offset by 1500kms^^with 
respect to the AGN. The y-axis has been reversed to facil- 
itate comparison with Figure El Object a is at the origin 
of the y-axis. 
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Fig. 5. Optical and near-IR spectra of object a, extracted 
with 1" X 1" and 0'.'68 x 1'.'2 apertures, respectively. 
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Line 






10'" X Pint 


/\ T'^T.^.TTTV.T 
— ■ r" vv riivi 


purest 




A 




erg s~^ cm~^ 


km s~^ 


A 


Lya 


4968.3±0.7 


3.086 


4.4±0.4 


1500±200 


120±20 


C IV A 1549 


6326.3±3.0 


3.084 


0.44±0.06 


800±300 


14± 3 


Hell A 1640 


6703.3±2.0 


3.087 


0.56±0.07 


700±225 


20± 4 


CIII] A 1909 


7791.6±3.5 


3.082 


0.47±0.07 


700±300 


24± 8 


[0 II] A 3727 


15280±100 


3.10 


3.9±0.3 




250±150 


H/3 


19840±100 


3.08 


0.7±0.2 




27±15 


[0 III] A 5007 


20370±100 


3.07 


11.5±1.0 




370±80 



° Deconvolved with the instrumental resolution. 



Table 3. Radio components of B3 J2330+3927. The 
8.4 GHz fluxes are normalized the northern component. 



Component 


RA(J2000) 


DEC(J2000) 


Normalized flux 


NVSS 


23''30™24!91 


39°27'11''4 




North 


23'' 30"* 24=836 


39°27'12'.'31 


1.00 


Centre 


23''30'"24!897 


39°27'11'.'54 


0.43 


South 


23''30'"24!922 


39°27'10''67 


0.23 



3.2. Where is the AGN? 

Our spectroscopic results identify object a as the most 
likely position of the AGN. This appears contradictory 
to the 8.4 GHz radio morphology shown in Figure |21 
and summarized in Table |3| The most straightforward 
interpretation of the VLA map is a core s urrounded by 
2 hotspots, as seen in many HzRGs (e.g. ICarilli et al.l 
llQQTtlPentericci et al.ll200o(i . The central component ap- 
pears marginally resolved with our 0'.'35 x 0'.'24 resolution, 
which looks lik e a carbon-copy of the z = 4.25 radio galaxy 
8C 1435-^635 l|Lacv et al.Lll994|) . The radio morphology 
therefore suggests that the central component marks the 
position of the AGN. Unfortunately, we do not have spec- 
tral index information on the individual radio components 
to examine this interpretation in more detail. 

Given the 0'.'4 astrometric uncertainty, we can state 
that the central radio component does not coincide with 
object a. If the nominal astrometry in Fig. [5] is correct, 
then the northern and southern radio lobes coincide with 
objects a and b, respectively. We can also exclude the pres- 
ence of a significant radio component north of object a be- 
cause the positions of the low-resolution radio data from 
the NVSS and from the VLA maps of Vigotti et al. (1999) 
(cross in Figure|2) are within 0'.'2 of the central radio com- 
ponent in our 8.4 GHz map. 

We propose two possible interpretations of the ob- 
served radio/near-IR configuration. First, the northern 
component may be the radio core of an extremely 
asymmetrical source, or a one-sided i et, much like the 
z = 4. 1 radio galaxy TN J1338-1942 JPe Breuck et al.1 
' ' hi' Il99^ have ar- 



ll999|) . lMcCarthv. van Breugel. fc K apah? fiyyilj have ar 
gued that the extended optical emission lines are bright- 
est on the side of the closer of the two radio lobes, which 
would be consistent with our spectroscopy (Fig. 01). The 
slight mis-alignment of the 3 components may then indi- 



cate a precession o f the radio jet, which i s often seen in 
FR II sources {e.g. iDennett-Thorpe et al.ll2002tl . 

The alternative interpretation is that the AGN is lo- 
cated at the position of the central radio component, but 
is eith er heavily obscured by dust {e.g. iReuland et all 
2003al). or not located at the centre of the host galaxy {e.g. 



Quirrenbach et all l200l|) . A potential problem with this 



interpretation could be the ionization of gas at object a, 
10 kpc away from the AGN. However, spatially extended 
CIV A 1549 and H ell A 1640 emission has been seen in 
sever al HzRGs {e.g. lOverzier et all liooH iMaxfield et all 

iooi). 



3.3. Possible synchrotron contribution at mm 
wavelengths 

Because B3 J2330-f-3927 is a powerful radio galaxy, it 
is important to estimate the possible contribution from 
synchrotron emission in the observed mm and sub-mm 
bands. We used the multi-frequency data of IVigotti et alJ 
(,1999) to determine the radio to submm spectral en- 
ergy distribution. We also ad ded the 365 MHz flux den- 
sity from the Texas survey l|Douglas et all h99^, our 
250 GHz MAMBO dete ction (§2.5), the SCUBA 353 GHz 
(850 Aim) detection of iReuland et all l|2003c(l . and the 
3mm (112.7 GHz) continuum emission seen in our PdBI 
data (see next section). We also include our tentative 
1.25 mm (230.5 GHz) detection, although we suspect that 
this emission has been over-resolved (see §3.5). 

Figure |^ shows this combined spectral energy distri- 
bution (SED). It is obvious that a synchrotron contri- 
bution would be negligible at 353 GHz (850 /im), but 
may contribute to the 113 GHz and 230 GHz detec- 
tions. A constant power-law extrapolation of the syn- 
chrotron spectrum from the 10.6 GHz data would predict 
a synchrotron flux density of 1.3 mJy at 113 GHz and 
0.6 mJy at 230.5 GHz. However, the synchrotron spec- 
tra of po werful radio gala xies often steepen at h i gh fre - 
quencies ijArchibald et all 12001: Andreani et all kOO'j] 
with spectral steepenings up to Aa ~ 0.5 l)Miirgia et all 
11999'). This suggests the constant power-law extrapola- 
tion should be considered as an upper limit. To esti- 
mate the 112.7 GHz contribution more accurately, we have 
fit the radio SEP, using the Continuous Injection model 
of IPacholczvkl l)l97fli) and the Single Injection model of 
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IJaffe k Perolal l)l973|) . Even with the latter model (which 
leads to an exponential steepening of the spectrum at high 
frequency), we found that the minimum flux density pre- 
dicted for the 112.7 GHz point is ~ 0.6 - 0.8 mJy. More 
observations between 10.6 GHz and 113 GHz would be 
needed to predict this contribution more accurately. 

Note that these values refer to the integrated flux den- 
sity of the radio source. However, Figure [51 shows that the 
radio source is resolved into three components, which ap- 
pear to have roughly similar flux densities at 8.4 GHz. 
Unfortunately, we do not have spectral index informa- 
tion of the individual components, but in most radio 
gala xies, the outer lobes have steeper spectra than the 
core ijCarilh et al.Lll997tlPmtericci et allboonl) . although 
sources with relatively steep spectrum cores have also been 
seen at high redshift ijAthreva et al.l Il997|) . The orienta- 
tion of the 112.7 GHz PdBI synthesized beam is roughly 
perpendicular to the radio axis, which means we should 
slightly resolve the synchrotron component at 112.7 GHz. 
Depending on the spectral index of the northern hotspot, 
this further reduces the synchrotron contribution at the 
position of the AGN to < 0.5 mJy. 



Observed Wavelength [m] 
1 0.1 0.01 0.001 

1000 — — — — n 




0.1 1 10 100 1000 



Observed Frequency [GHz] 

Fig. 6. Radio to submm spectral energy distribution of 
B3 J2330+3927. The sohd hue is a modified black body 
curve with T = 50 K and /3 = 1.5, normalized at 
352 GHz. The 230.538 GHz PdBI measurement severely 
over-resolves the thermal dust emission. 



3.4. Identification of the CO emission 

Figure 13 shows the 3mm PdBI data, integrated between 
112.515 GHz and 112.723 GHz, which represents the fuU 
width of the redshifted CO J = 4 — 3 line (see below). The 
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Fig. 7. Integrated CO J 4 — 3 emission towards 
B3 J2330-1-3927 overlaid on the Keck/NIRC X-band im- 
age. Contour levels start at 2 a, and increase by 1 cr, with 
(T = 0.2 mJy Beam~^ (negative contours are dotted). The 
3 crosses mark the positions of the radio components in 
Fig. 121 The FWHM of the synthesized beam is shown in 
the bottom left corner. 

emission clearly peaks at the position of object a, which 
is the most likely position of the AGN (see § 3.2). The 
emission does not follow the extension of the radio source, 
confirming that a possible synchrotron contribution would 
be minimal. 

Figure |H1 shows the PdBI CO J = 4 - 3 spectrum of 
B3 J23304-3927 extracted at the position of object a. We 
show the unbinned spectrum, and a spectrum binned to 
105 kms^^. A broad emission line is clearly seen around 
112.6 GHz, with a weak underlying continuum emission. 
We fit a Gaussian profile with an underlying continuum 
emission to determine the line and continuum parameters. 
Table 0] gives the fit parameters for 5 different bin sizes 
used. As long as the bin size does not exceed 40 MHz, 
the derived parameters are not very sensitive to the bin 
widths. 

In the previous section, we estimated the contribution 
from synchrotron emission to be < 0.5 mJy, so the broad, 
~2 mJy line appears securely detected. We interpret this 
emission line as CO J = 4 — 3atz = 3.094, i.e. offset by 
~500 kms^^ to the red from the optical redshift, deter- 
mined from the Hell A 1640 line. Such offsets have also 
been seen in high S /N CO iden tification s of high redshift 
quasars {e.g. Guillotei^TT^n i]. J999; Cox et all Efio2^ . 
The ^0.3 mJy continuum emission seen in the 1.3 mm 
spectrum (Fig. ^ is consistent with the extrapolation of 
the thermal dust spectrum (Fig. (SJ, but it is also likely 
that there is a synchrotron contribution from the north- 
ern radio component, especially if the latter would have a 
flatter spectral index than the integrated radio emission. 
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The total velocity width of the CO-line is 
~500knis^^(sce Fig. (S)), and the integrated 
flux over these five 106.5 kms~^ channels is 
5*00^^ = 1-3 i 0.3 Jy kms"^. Note that this value has 
been corrected for the underlying continuum emission 
0.3 mjy). This represents one of the faintest un-lensed 
detections of CO emission to date. 

To determine if the CO emission in spatially resolved, 
we have subtracted different source models (point source, 
circular or elliptical Gaussians) from the UV-data. We 
find no clear evidence that the integrated CO emission 
is spatially resolved, despite the suggestion in Fig [3 that 
the emission appears resolved towards object c. The latter 
may be due to the low S/N of the outer contours. Note 
that due to our limited frequency coverage, we may have 
missed CO emission from components with a large velocity 
offset. For example, object b is known to have a velocity 
offset of 1500 kms~^(Fig. so any CO emission from 
this object would not have been detected in our observa- 
tions. Deeper integrations with a wider frequency coverage 
would be needed to examine the spatial structure of the 
CO emission. 



Velocity Offset (z = 3.087) [km/s] 
-500 500 1000 




Observed Frequency [GHz] 



Fig. 8. CO J = 4 - 3 spectrum of B3 J2330-I-3927 ex- 
tracted at object a. The top panel shows the unbinned 
spectrum, while the bottom panel has been re-binned with 
Ai/chan = 40 MHz (106.5 km s-i). The velocity scale is 
with respect to the central frequency ofCOJ = 4 — 3at 
z = 3.087. The solid line is a Gaussian fit with an under- 
lying continuum emission. 



Table 4. Gaussian fitting coefficients of the CO line in 
object a 



Bin width 


# bins 


'S'cont 


'S'pcak 


^centre 


AVpwHM 


MHz 




mJy 


mJy 


GHz 


km/s 


2.5 


270 


0.36 


2.51 


112.620 


207 


10.0 


67 


0.33 


2.55 


112.620 


211 


20.0 


33 


0.36 


2.47 


112.623 


208 


40.0 


16 


0.28 


2.51 


112.621 


223 


80.0 


8 


0.52 


2.22 


112.614 


208 



3.5. Identification of the spatially extended thermal 
dust emission 

From the radio to submm spectral energy distribution 
(Fig. EJ, it is obvious that the interferometer data under- 
estimate the dust emission at 230 GHz. We now examine 
this possibility in detail using our 230.5 GHz PdBI map 
(Fig.©. 

Extrapolating the 353 GHz (850 fim) SCUBA detec- 
tion (14.15 ± 1.65 mJy) using a blackbody with T = 
50 K predicts an integrated dust emission fiux density of 
4.2 mJy at 230 GHz and 0.45 mJy at 113 GHz. Figure El 
shows the central part of our 230.5 GHz PdBI map. We 
find a marginal detection of dust emission with an in- 
tegrated fiux density of 1.2±0.9 mJy. This appears in- 
consistent with the SCUBA detection. Such a discrep- 
ancy has also been seen in high redshift quasars, where 
the PdB interferometer data underestimate the flux by 
as much as an order of magnitude when compared to 
IRAM 30m single dish rneasurements at the sam e fre- 
quency tauillotea u et all Il999t lOmont et all. l200l[l . The 
most likely explanation is that the dust emission is ex- 
tended with respect to the much smaller interferometer 
synthesized beam. 

Although our 230 GHz detection is marginal at most, 
it is concentrated near the peak of the CO emission and 
the most likely position of the AGN (see §3.2), support- 
ing the reality of this detection. We conclude that there 
is evidence for a two-component spatial distribution of 
the dust emission: flrst, a compact component, which may 
be heated by the AGN, and second, a component signifi- 
cantly larger than the 1'.'9 x 1'.'5 PdBI synthesized beam 
at 230 GHz. 

The fact that the ^0.3 mJy continuum emission seen 
in the 113 GHz spectrum (Fig. ^ is close to the extra- 
polation of the thermal dust emission (despite a possible 
synchrotron contribution) suggests that the 3'.'2 x 2'.'3 syn- 
thesized beam does not over-resolve this dust emission. In 
theory, we could determine the scale size of the thermal 
dust emission by combining the zero-spacing JVIAJVIBO 
flux density (scaled to 230 GHz) with the PdB interferom- 
eter data fFig. 1101) . However, due to the low signal/noise 
of our data, we can only constrain the scale size between 
0'.'5 and 5". 
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Fig. 9. IRAM PdBI 230.5 GHz (1.3 mm) map of 

B3 J2330+3927 overlaid on the Keck/NIRC i^-band im- 
age. Contour levels start at 1 a, and increase by 1 a, with 
(T = 0.45 mJy Beam^^ (negative contours are dotted). 
The 3 crosses mark the positions of the radio components 
in Fig. 13 The FWHM of the synthesized beam is shown 
in the bottom left corner. 



100 200 
Baseline Radius [m] 

Fig. 10. Visibility amplitude versus baseline radius in the 
the uv-plane for the 230 GHz PdB interferometer data 
(black/thin points), and the MAMBO/IRAM single dish 
data (thick/red point) interpolated to 230 GHz. The solid 
line represents a point source of 1.8 mJy, the dotted/blue 
and dashed/magenta lines a circular Gaussian source with 
a FWHM of 0'.'5 and 5", respectively 



3.6. No HI 21 cm absorption 

Despite our deep and repeated WSRT observations, 
we have not been able to convincingly detect HI 
21 cm absorption against the radio source. This is 
consistent with unsuccessful WSRT search es in 6 
HzRGs (jRotteering. de Bruvn. fc Pentericcl Il999h . 
To date, HI 21 cm absorption has been solidly 
detected in only one HzRG, viz. B2 09 02+34 
iz = 3.3909 llTson. Bagri. CornwellL l1991^ 
iBriggs. Sorar. fc Taramopoulost 1 19931) . This raises 
the question as to why B2 0902+34 stands out. One clue 
may c ome from its bizarre radio morph o logy llCarillil 
Il995|) . iRottgering. de Bruvn. fc Pentericcil l)l999l) argue 
that the HI is located in a kpc-sized disc, which can 
be seen in absorption due to a fortuitous orientation 
of the radio source. They further suggest that a high 
core fraction may characterize such orientations: 15% 
for B2 0902+34, compared to 0.6% and 2.9% for the 6 
HzRGs where no 21 cm absorption has been detected. 
However, in B3 J2330+3927, we find a core fraction of 
60% or 26% (see Table OJ, depending if we adopt the 
northern or central radio component as the core (see 
§3.2). We conclude that the present sample of 8 HzRGs 
where HI 21 cm absorption has been searched is clearly 
too small to address the origin of the absorbing neutral 
gas. 



4. Mass estimates 

The detection of the different emission and absorption 
components in B3 J2330+3927 suggests it is a massive 
system. In this section, we use the above data to derive 5 
independent mass estimates of the gas and dust compo- 
nents. We discuss each component separately, and sum- 
marize our results in Table El We note that several high 
redshift CO detections have been helped by gravitational 
lensing (for reviews, see 'Combes. M aoh. fc Omoii^. Il999t 
|£!ox, Omont fc Bertoldi, 2002b). However, from our radio 
and optical/near-IR imaging, we find no indications that 
B3 J2330+3927 could be amplified by gravitational lous- 
ing, and we do not consider this in the following. 

4.1. Molecular gas mass 

We can estimate the molecular gas content from our CO 
line detection by making an assumption for the standard 
conversion factor Xqo = M {H2) / L'qq{1 — 0), where 

= (l^)'^coA^^o-b.^i(l + (1) 

is the CO line luminosity me as ured i n K kms~^ pc^ 
llSolomon. Downes. fc Radfordl ll992|) . The stan- 
dard value for Galactic molecular clouds is 
A^CO = 4M0(Kkms^^pc^)^-^. However, there are strong 
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indications from studies of extreme starbursts in ultra lu- 
minous infra-red galaxies (ULIRGs) that this value should 
be considerably smaller for dense starbursts and AGN. We 
shall theref ore assume Xgn ^ 0.8 M0(Kkms~^pc^)~-'^, 
as given bv lDownes &: Solomon (1998). 

Because we have observed the CO J = 4 — 3 transi- 
tion we also need to esti i nate t he (4 — 3)/(l — 0) line ratio 
r43 . IPapadopoulos et all l)2000|) have used a large velocity 
gradient code to predict the r43 ratio in two radio galaxies 
at similar redshifts as B3 J2330+3927. Using an interstel- 
lar matter environment similar to that of a local average 
starburst, the y predict = 0.45. 

Following ISolomoiK Downes. fc RadfordI l|l992|) . we 
find i^o(4 - 3) = 3.9 x lO^^h^^ K kms'^ pc^. Inserting 
the values given above, we find 



M(H2) 



^co 

f34 



This value is slightly lower than the ones found in the pre- 
vious two HzRGs, 4C 60.07 and 6C J1908-h7220 (2.4 and 
1.35 X IO^^Mq, respectively, in our adopted cosmology). 

4.2. Dynamical mass from CO 

Our PdBI CO data does not have sufficient S/N to de- 
rive detailed spatial and spectral information needed to 
fit models of rotating disks, such as those often seen in 
local ULIRG s and recently also at high redshift (z = 2.8; 
iGenzel et al.L E(103i) . However, we can derive a rough up- 
per limit on the dynamical mass by assuming that the 
size of a possible rotating disk is smaller than the resolu- 
tion of our PdBI 113 GHz observations. Provided that the 
CO is effectively located in such a disk, this is probably 
close to the actual extent, given our estimated scale size of 
the dust emission, which is often correlated with the CO 

emission. 

Following IPapadopoulos et al.l l)2000|) . the dynamical 
mass Mdyn is then given by 



M, 



^ dyn 



1.16x10^ 



with (Sobs the observed flux density at a wavelength where 
the SED is dominated by thermal dust emission, Kci(t'rest) 
the rest- frequency mass absorption coefficient, B(v^cst, ?d) 
the Planck function for isothermal emission from dust 
grains radiating at a temperature Td, and the lumi- 
nosity distance. 

We need to make some assumptions for and 
Kd(t'icst)- Because we observe the steep Rayleigh- Jeans tail 
of the dust distribution, a small change in the adopted 
dust temperature may cau se a significant change in the 
derived dust mass (see e.g. iHughes. Dunlop. fc Rawlingsl 
Il997t iMcMahon et all \l99^ . The uncertainties in both 
the SCUBA 850 and the MAMBO 1.3 mm flux den- 
sities cannot constrain Td tightly (see Fig. EJ. We will 
therefore adopt Td = 50 K, w hich has been used for most 
other high redshift AGN (e.g. IPapadopoulos et al.Ll200(i 



FWHM 



100 km s 



-1 



A (3i,^)-Mo,(2) 



where i is the inclination of the disk. Adopting T < 25 kpc 
and AVfwhm = 220 kms^^, we find 

Mdyn < 6.4 X 10i°(sin2i)^iMQ. 

Although our estimates of M(H2) and Mdyn are very un- 
certain, this suggests that a substantial fraction of the 
total mass is in the form of molecular gas. 



4.3. Dust mass 

FigureHshows that the 850 fim flux of B3 J2330+3927 is 
almost completely dominated by thermal dust emission. 
We can therefore estimate the dust mass Aid using: 



Md = 



(1 + z)Kd(t'rost)S(i^rcst,Td) 



(3) 



lArchibald eral].l200ltlbmont et al.U200l|) . and OTOvides ; 
good fit to the thermal part of the SED of B3 J2330-t-3927. 

The rest- frame mass absorption coefficient «;d(t'rost) 
can be found by extrapolation from other fre quencies us- 



i Chini. Kriigel. fc Krevsal Il98rf>. This emis- 



mg Kd (X 

sivity index (3 can be determined from multi-frequency 
observations covering the entire dust emission spectrum. 
However, such observations are scarce for high redshift 
objects, and assume that the dust emission is dominated 
by a single temperature component. Because both (3 and 
Td determine the shap e of the dust emission spec trum, 
their values are coupled llPriddev fc McMahoniEffll . The 
publis hed values of var y between (3 = 1.3 for Td = 
36 K (Dunne et all 120 00") and /3 = 1.95 for Td = 41 K 
(jPriddev fc McMahonL i2001.1 . We shall a dopt the inter- 
medi ate value of /3 = 1.5 for Td = 50 K ijBenford et all 
[1999). These values are consistent with our 113 GHz to 
352 GHz measurements (Fig. 01. 

The value of Kd depends on the composition of the 
dust grains. The estimates of Kd normalized to 850 /im 
using (3 = 1.5 vary b etween Kd(850/im) = 0.04 m^kg^^ 
Draine fc Lei Il984l) an d Kd(850/^m) = 0.3 m^kg-^ 



Mathis fc WhiffenLll989l) . We shall normalize to an inter- 
mediate value of «;d(850/Ltm) = 0.11 m^kg"^ (|HildebrandL 
Il983l) . 

Using these values and S'sso^un 
(|Reuland et alil2003cj) . we find 



14.15 ± 1.65 mJy 



Md 



9 X lO^/igt 



0.11 



T(Td) M, 



0j 



Kd(850/im)^ 

with T(Td) = (exp[/ii//fcTd] - l)/(exp[/ii//fc • 50K] - 1). 

Comparing with the molecular gas mass estimate, this 
yields a gas-to-dust ratio 

M(H2) _ 850 / 0.11 \ 
Md ~ T(Td) Ud(850Mm)j ' 

This is comparable with values found in other 
HzRGs fPapadopoulos et al.', '2000*) and ULIRGs 
ijSanders. Scoville. fc Soifer , _1991). Alternatively, we 
can also compare the CO and 850/im fiuxcs, yielding a 
comparison that is independent of the assumed param- 
eters. We find Ay5(CO)/S'800Mm = 92 kms-\ again 
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si milar to the values found in the 2 HzRGs observed 
bv IPapadoDOulos et alJ ll2000l) the 3g hmit in 53W002 
ijAU oin. Barvainis^_fc_Gjjjlloteau|^ 200Cil. a nd high redshift 
quasars (^Cox. Omont fc Bertoldi . l2002b(l . This indicates 
a standard gas/dust ratio in B3 J2330+3927. 

4.4. Ionized gas and associated H I absorption in Lya 

Figure 21 shows the complex velocity structure of the Lya 
line. The emission spatially extends well beyond the ra- 
dio l obes. This has been seen in several other HzRG s 
(e.g. iDe Rreuck et all Exvil IVillar-Ma.rtfn et a].L Einj. 
and excludes shocks as the ionization mechanism at this 
position. Because we do not detect a discontinuity in the 
Lya profile at the position of the radio lobes, this implies 
that photo-ionization by the central AGN is the dominat- 
ing ionization mechanism of the Lya line. 

Velocity Offset (z = 3.087) [km/s] 
-2000 -1000 1000 2000 



°< 20 




4940 4960 4980 

Observed Wavelength [A] 



5000 



Fig. 11. Lya velocity profile of B3 J2330+3929 
(solid/black histogram). The solid/blue line is the model 
consisting of a Gaussian emission (dashed/red line) 
and two Voigt absorption profiles with the indicated 
parameters. The dotted/green histogram is the scaled 
He II A 1640 profile. The horizontal bar indicates the full 
velocity width of the CO 4 - 3 line (Fig.|Hl). 

Figure [TTl zooms in on the Lya line in the Keck/LRIS 
spectrum (extraction width 1" x 1"). To check if the ve- 
locity profile can be explained by pure velocity structure 
without absorption, we have attempted to fit the Lya 
profile with two or three Gaussian components, but this 
does not provide a good fit between the two peaks, or on 
the blue wing. An alternative interpretation is that the 
double-peaked structure is caused by associated HI ab- 
sorption. Such absorbers are often seen in the Lya lines 



of HzRGs ( e.g. Ivaii Oiik et all ITQQTat iDe Breuck et all 
l2nonbH,Tarvis et SilLhOQ^ . 

We therefore modeled the profile with a Gaussian emis- 
sion profile, supplemented by two Voigt absorption com- 
ponents (a second component is needed to fit the non- 
Gaussian blue wing). Figure [TTl presents our best fit, ob- 
tained through a minimization. We find that the profile 
can be closely represented with this model. Although the 
central wavelength of the Gaussian emission was a free 
parameter in our fit, it corresponds very well with the 
redshift obtained from the He II A 1640 line, providing a 
strong physical base for our model. 

We can estimate the mass of the associated HI ab- 
sorber using 



MiRl) = mpRl^,N kg, 



(4) 



where nip is the proton mass, i?abs is the size of the ab- 
sorption system, and N is the column density derived from 
the Voigt profile fitting. In convenient astrophysical units 
ijvan Oiik et alJ .[l997a). this becomes: 



Af(HI) 



M, 



0j 



(5) 



where R35 is in units of 35 kpc, and iVig is in units of 
10-i9cm2. 

Because the red absorber has a ^10000 times higher 
column density {N — 8.7 x 10^* cm^^) than the blue one, 
we shall only consider the red absorber in the following. 
From Figure ^ we obtain a lower limit of > 2" for the 
full spatial extent of the red absorber. Assuming conser- 
vatively a size of i?abs ~ l^/i^^^ kpc, we derive a total mass 
estimate of the absorber 

M(HI) > 2 X IO^Hq^Mq. 

We can also use the Lya emission to estimate the mass 
of ion ized hydrogen M(HII). Following McCar thv et alJ 
l|l99C|) , we assume pure case B recombination at a temper- 
ature of T = lO'' K. Using the Lya emission corrected for 
the associated HI absorption (dashed line in Figure [TT]l . 
we obtain Af(HII) using 



Af(HII) = lO9(/_5i44l^7o)'/'M0, 



(6) 



where /_5 is the filling factor in units of 10~^, L44 is the 
Lya luminosity in units of 10'*'' ergs s~^, and V70 is the 
total vol ume in units of 10^° cm '^. Assuming a filling factor 
of 10-5 (|McCarthv et al.L ll99nV and a volume of 6" x 2" x 
6" = 1.4 X 10^9 /igg^ cm^ (see Fig. El, we find 



2.5 X lO^hQ^^^Mr. 



A/(HII) ~ Z.U X lU Jlgg iVlQ. 

This valu e is consis tent with those found in other HzRGs 
{e.g. van Oiik et al... 1997a). 

4.5. HI 21 cm absorption against the radio source 

We can determine an upper limit of the H I mass from our 
non-detection of the HI 21cm line. The optical depth is 
defined as 



r = - ln(l-A^/(5cf)), 



(7) 
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Table 5. Mass estimates of gas and dust components in 
B3 J2330+3927, subject to uncertainties discussed in §4. 



Component Mass estimate Determined from 

Me 

7 X 10"" CO line 

Dynamical < 6.4 x 10"'(sin2 i)"" CO line 

Dust 9 X 10*^ 850/im continuum 

HI 2 X 10*^ Lya absorption 

HII 2.5 X 10* Lya emission 

HI < 2 X 10** HI 21 cm absorption 



where AS is the H I absorption flux, S the continuum flux 
and Cf covering factor that is assumed to be 1. Assuming 
a 3ct upper limit, we find r < 1.3%. 

A hmit to the integrated colu mn density ca n then be 
derived from the expression {e.g. lRoh1ffJ.ll9s4 : 

Nil ~ 1.83 X 10^^T,pi„y" rdv cm'^, (8) 

where Tspin is the spin temperature in Kelvin, and we 
assume a profile width of 300km s~^. HI studies of 
high-redshift damped Lya s ystems (jCarilli et all Il996t 
iKanekar fc Chengaluil l2003|) have found values as high 
as Tgpin = 1000 K, but there also appears to be a trend 
for lower Tspin in the most luminous galaxies. We therefore 
adop t the Tsnin = 100 K value found in Galactic clouds 
{e.g. iBranm Walt^Th^lTflfli. This yields 

iVn < 9.3 X 102°(rspin/100) cm^^. 

Assuming that the absorption occurs only against the ra- 
dio core, we obtain an estimate of the scale-length of the 
absorber of R^hs ^ 0'.'5 « 4 kpc. Using equation (4), this 
yields 

M(HI) < 1 X lO^/lgg^M©. 

5. Discussion 

5.1. The interaction between the different gas 
components 

We now examine the relationship between the different 
components summarized in table [3 The similarity of ob- 
jects b and c in both size and iiT— band magnitude suggests 
that object c may also be part of the same system. We 
may therefore witness the accretion of the smaller com- 
panion object b (and possibly also object c) by a mas- 
sive galaxy hosting an AGN (object a, but see §3.2), sur- 
rounded by a large reservoir of molecular gas and dust. 
Such gas reservoirs are often directly seen in the huge 
Lya haloes , which extend out t o >100 kpc from the 
AGN ( e.g. Ivan Oiik et all Il996l: IVenemans et~all l2002t 
IVillar-Martm et all l20oTlReuland et all l2003l| ^ From 
Figure ^ we find that the extent of the Lya halo in 
B3 J2330+3927 is --6" (50 kpc), but our spectrum is not 
very deep, so this value should be considered as a lower 
limit. 



The peaks of the CO, dust and Lya emission are all at 
object a, suggesting that the ionized hydrogen and the CO 
and dust emission are co-spatial. Figure^Jshows that the 
CO emission also coincides very closely in velocity space 
with the main associated H I absorber. Because these com- 
ponents are both co-spatial (at object a), and have similar 
velocity shifts, this strongly suggests that we detect CO 
emission from the same cloud that absorbs the Lya emis- 
sion. Our mass estimate from the associated Lya absorber 
is ~ 3500 times lower than the one obtained from the CO 
emission (Tabled . Although both estimates are uncertain 
to almost an order of magnitude, this difference does indi- 
cate that the associated Lya absorber probes only a tiny 
part of the total gas mass, either because most of the hy- 
drogen is in molecular form, or due to the requirement of a 
strong background source to detect the H I in absorption. 

Our multi- wavelength data of B3 J2330-I-3926 show for 
the first time a direct relationship between the associated 
HI absorption and CO emission in a HzRG. High reso- 
lution spect roscopy of the Lya and CIV A 1549 li nes of 
two HzRGs teinette et all l2000t IjaTvis et al.Ll2002^ have 
shown that the absorbers cannot be co-spatial with the 
emitting gas. The absorption should therefore occur fur- 
ther outward, possibly in a shell pressur ized by the ex- 
panding radio source ijBinette et all EoOQl . Our detec- 
tion of CO emission from the absorbers suggests that 
this cloud may be quite extended, but it cannot be pri- 
mordial, as it is already chemically e nriched. This is als o 
consistent with the recent findings of lBaker' et all ll2002l) . 
who found a strong correlation between the equivalent 
width of CIV A 1549 absorbers and optical continuum 
slope in radio-loud quasars. They interpret this relation 
as evidence for the presence of dust in the associated 
CIV A 1549 absorbers, which reddens the optical con- 
tinuu m. The gas and dust emissi on often trace each other 
{e.g. IPaoadoDOulos et an i200rf). suggesting that the as- 
sociated absorbers seen in both quasars and radio galaxies 
probe the large gas and dust reservoirs surrounding AGNs. 

The presence of CO from these reservoirs implies 
that they contain processed material. This metal emis- 
sion is sometimes also seen in emission in extended 
haloes of CIV A 1549, Hell A 1640 and NV A 1240 
dMaxfield et all l2002t IVillar-Martm et all l2003l) . They 
may have been deposited by previous merger events, 
or by star burst-driven superwinds at even higher red- 
shifts (e.q iTaniguchi et all l2001at IPawson et all 120021: 
Aiiki et all l2002t iFurlanetto fc Loebl l2003l) ~ndeed. 
Taniguchi. Shiova. fc KakaziJ l|2001m discuss evidence 
that a superwind occurs in a z = 3.09 "blob" of Lya 
emission which a, l so eni its strongly at submm wavelengths 

a oman et alll200lh . Such objects may be related to 
HzRGs, e.g. when observed during a radio-quiet phase. 

To summarize, our observations of B3 J2330-I-3927 
provide further evidence that HzRGs are mas- 
sive galaxies, located in a dense, interacting en- 
vironment. Deep near- IR imaging wit h Kec k 
(van Brcugcl et al.', 'l998^ IPe Breuck et all l2002l) . 
HST (Pnntericci et al., .20011 . and ground-based adaptive 
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optics l|Steinbring. Crampton. fc Hutching^ l2002j) shows 
that the host galaxies are often surrounded by fainter 
clumps, hke objects b and c in B3 J2330+3926. These 
clumps may well be concentrations within the large 
gas/dust reservoir surrounding the AGN, which will 
eventually merge with the massive central galaxy hosting 
the AGN. 

5.2. Star formation rate 

Having found evidence that the host galaxy of 
B3 J2330+3927 is surrounded by a large gas reservoir, 
which has already been chemically enriched, it is of in- 
terest to estimate the ongoing star formation rate in this 
galaxy. We cannot use the rest-frame UV and optical emis- 
sion lines, because they are mainly ionized by the AGN, 
and not by massive stars. Similarly, the UV continuum 
emission is likely t o be contaminated by a scattered quasar 
component {e.g. IVernet et al.l l2001a|) . We therefore de- 
termine an upper limit to the global star formation rate 
of the entire system using the far-IR (FIR) dust emission. 

We can determine the total FIR luminosity by inte- 
grating the thermal spectrum: 

poo 

LFiR = 47rAfd/ Kd{'^)B{iy,Td)diy, (9) 
Jo 

yielding: 

^FiK = ^ ■ ''^np ^mp+ 4]M„ (10) 

where T and C are the Gamma and Riemann C functions, 
respectively. 

Substituting equation (3) and x = hv-^cst I kT^, we find: 
Lfir = 47rr[/3 + 4]C[/3 + A]DIx~^P+'^^ (e^ - l)S'obsJ^obs, (H) 
For the values adopted in this paper, we find: 
Lfir » 3.3 X lOi^/ig/Lg. 

In §3.5, we have argued that the dust emission prob- 
ably consists of a central component near the AGN, and 
a spatially extended component. This would imply that 
the power source of the FIR emission is a combination of 
direct heating by the AGN (heating the central compo- 
nent) and by recently formed massive stars in a starburst 
(heating the extended component). It is unlikely that the 
AGN can power the dust emission out to several tens of 
kpc, while the detection of several companion objects in 
the optical and near-IR images suggests the presence of 
stars out to ~40 kpc. 

Following lOmont et all (|20o3), we can then calculate 
the star formation rate 



SFR = 5MF<5sB(Lf /fl/lO^^Lo) Moyr-\ 



(12) 



where 5mf is a function of the stellar mass function, and 
^SB is the fraction of the FIR emission heated by the star- 
burst. Assuming conservatively 5mf = 1, we find 

5Fi? « 3000(5sb/^65 Moyr-i. 



IPapadopoulos et all l)200Cl|) report star-formation rates 
up to 1500 M0yr~^ in the 2 other HzRGs where CO emis- 
sion has been detected. Because ^sb is likely to be sig- 
nificantly smaller than 1, our values are probably lower. 
Nevertheless, they are still extremely high, indicating that 
we are witnessing a major starburst phase, possibly trig- 
gered by the interaction and/or merging of the different 
nearby companion objects seen in the if— band image. 

6. Conclusions 

The main results from our multi-frequency observations 
of B3 J2330-I-3927 are: 

• We detect the CO J = 4 — 3 line with a full ve- 
locity width of ^^500 kms~^ at the position of the AGN 
host galaxy. The line is centered at z = 3.094, much closer 
to the central velocity of the associated HI absorber in 
Lya than to the systemic redshift determined from the 
Hell A 1640 line. This strongly suggests the CO emis- 
sion and associated HI absorption both originate from 
the same gas reservoir surrounding the HzRG. 

• We find a substantial discrepancy between our 
250 GHz single dish and 230 GHz interferometer measure- 
ments of the thermal dust emission. This strongly suggests 
that the dust emission is distributed at spatial scales larger 
than the 1'.'9 x 1'.'5 synthesized interferometer beam. The 
combined uv-plane data can be fit by a circular Gaussian 
profile with a size between 0'.'5 and 5". 

• The optical/near-IR spectroscopy suggests that the 
AGN is located at the brightest if— band source (ob- 
ject a). This is supported by the CO and dust emission, 
which both peak at the same position. However, the radio 
morphology suggests that the AGN is located at the cen- 
tral radio component, 1" south of object a. If the AGN is 
located at object a, the radio source is extremely asym- 
metric. Alternatively, the AGN may be located at the cen- 
tral radio component, and could be either heavily obscured 
in the optical/near-IR, or could be located off-centre from 
the host galaxy. However, this alternative appears unlikely 
given the position of the dust/CO peak. 

• We obtain a t < 1.3% limit on the H I 21 cm absorp- 
tion against the radio source. This represents the seventh 
non-detection out of 8 radio galaxies where H I 21 cm ab- 
sorption has been searched. 

• We estimate the mass of the different components 
in B3 J2330-t-3927, which indicates that the CO emission 
traces the largest fraction of the total mass in the system. 

In summary, our multi-wavelength observations sug- 
gest the presence of a large gas and dust reservoir sur- 
rounding the host galaxy of B3 J2330-I-3927. For the first 
time, we have found strong indications that these differ- 
ent observations may trace the same material, both in 
emission and in absorption. The detection of CO emission 
therefore implies that a substantial part of this reservoir is 
not primordial, but must have been previously enriched, 
possibly by a starburst-driven superwind. 
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